Heat Conduction in Low Dimensions: Prom Fermi-Pasta-Ulam Chains to 

Single- Walled Nanotubes 
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Heat conduction in 1-dimensional anharmonic systems is anomalous in the sense that the con- 
ductivity k scales with a positive power of the system size, k ~ L a . In two dimensions, previous 
simulations and theoretical arguments gave a logarithmic divergence. For rectangular systems of 
size L« x L± there should be a cross-over from the 2-d to the 1-d behaviour as the aspect ratio 
r = Lu/L± increases from r = 1 to r > 1. When taking periodic boundary conditions in the trans- 
verse direction, this should be of direct relevance for the heat conduction in single-walled carbon 
nanotubes. In particular, one expects that k nanotubes of diameter R should conduct heat better 
than a single nanotube of the same length and of radius kR. We study this cross-over numerically 
by simulating the Fermi-Pasta-Ulam model. Apart from giving a precise estimate of the exponent 
a, our most intriguing results are that the divergence does not seem to be logarithmic in d = 2 but 
also power-like, and that the cross-over does not happen at a fixed aspect ratio. Instead, it happens 
at r — r* with r* — » oo for L — » oo. 
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After years of intense studies it is now clear that heat 
conduction in typical anharmonic systems is anomalous 
in low dimensions [fjj. In particular, in 1-d systems like 
the Fermi-Pasta-Ulam (FPU) /3-model the diatomic 
Toda chain § , or in 1-d hard-particle gases with alter- 
nating masses @, |], [||(| the heat conductivity k, defined 
by J(x) — k\7T(x), scales as 



(1) 



where L is the size of the system. Notice that we have 
assumed here that T(x) exists, i.e. that local thermal 
equilibrium (LTE) is established in the limit L —> oo, 
which is non-trivial but seems to be true for these 
systems. 

The main reason for this anomalous behaviour is the 
fact that the free path length of phonons with long wave 
lengths A diverges for A — > oo. Thus soft modes propa- 
gate nearly ballistically. In higher dimensions this is also 
true, but due to the enlarged phase space for other modes 
it is restricted to such a small region that k remains fi- 
nite. In 1 dimension there are just not enough other 
modes with which soft phonons could interact to make 
k finite. This argument requires of course that there is 
an acoustic phonon branch, i.e. that soft (Goldstone) 
modes exist because translation invariance is not broken. 
Thus this argument does not apply to charge conduction 
where the electrons move in an external potential braking 
translation invariance. 

For 2-d systems, simulations of FPU lattices of rather 
modest sizes || indicated a logarithmic divergence n ~ 
InL. This was also supported by mode coupling theory. 
Indeed, a logarithmic divergence should not be too sur- 
prising in view of the logarithmic divergence of transport 
coefficients in 2-d hydrodynamics || due to long time 
tails. In H it was also observed that the cross-over from 
the 1-d to the 2-d scaling in rectangular systems happens 
at surprisingly large aspect ratios. Consider a system of 
size Ln x L± where the length Lh is parallel to the heat 
flow. The aspect ratio is defined as r = Ln/L±. For the 



sizes studied in ||, the cross-over happened at r w 10, 
but no detailed study was made. 

The main aim of the present paper is to present de- 
tailed simulations of the FPU /3-model on much larger 
lattices than in ||] , in order to study in detail the asymp- 
totic behaviour of 2-d systems and of the cross-over. But 
before doing so, let us point out that this study is not 
entirely academic. It is of immediate experimental and 
maybe even of technological importance. 

The lattices which we shall use have periodic boundary 
conditions in the transverse direction. Therefore, we ac- 
tually study the heat conduction through (single-walled) 
tubes. By far the most important nanotubes are those 
made of carbon [jl(| . Their heat conduction is dominated 
by phonons. They can be made with very few lattice de- 
fects, such that their conductivity is basically controlled 
by phonon-phonon interactions as in the FPU model. In- 
deed, this conductivity has been both calculated |[l| and 
measured [^2| , and is found to be very large. It was thus 
suggested that carbon nanotubes should be the ideal ma- 
terial to carry away Joule heat in the next generation 
of integrated circuits which will use nanoscale structures 
and for which cooling will be a major problem. Since 
most molecular dynamics simulations and measurements 
were done for fixed tube lengths, no length dependence 
of k was seen in them. The only exception is Ref. |l3}| 
where a clear length dependence of k was seen for the 
narrowest tubes, but not for wider ones. 

In all simulations of carbon nanotubes the realistic 
Tersoff-Brenner potential |l4j was used. If, as we claim, 
the size dependence of n is universal, its main features 
(and in particular the exponent a in Eq. ([[])) should be 
independent of the potential. There will be of course un- 
known scale factors when transferring our results from 
FPU systems to carbon tubes, nevertheless we should 
be able to use our results directly for an experimentally 
accessible system. 

In our FPU simulations we used a square lattice with 
sites indexed by integer vectors n = (i,j) with 1 < i < 
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FIG. 1: Temperature profiles for lattices with N x = 128 and 
N y = 1,3, and 32. For these simulations we used a smaller 
step size dt — 0.025. 
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FIG. 2: Normalized kurtosis profiles from the same runs used 
also in Fig.l. 



N x and 1 < j < N y . We use periodic b.c. in j. The 
lattice constant, the particle mass, and the (3 parameter 
are all set to unity. Thus N x = Lu and N y = L± and, 
apart from the heat baths at i = 1 and i = N x , the 
hamiltonian is 



1 



H = - 



2 E^ 



<n,m> 



,) 2 + 



(2) 

Notice that identically the same form (with N y = 1) can 
be used also for 1-d chains. To measure a heat flux we 
keep the left boundary (i = 1) at a lower temperature 
T/i = 10, and the right boundary (i — N x ) at T c = 6. 
These are fairly large in order to obtain fast equilibra- 
tion (similar temperatures were used in ^|). These ther- 
mostats are implemented as Nose-Hoover |J heat baths 
with response times = 1. Thus the equations of motion 
for the particles in the boundaries are modified to 



dH 



Up* 



with 



P 2 

f — — - 1 



(3) 



Since our overall accuracy was mainly limited by statis- 
tical fluctuations and not by integration errors, we used 
mainly a simple leap-frog integrator [^) with a large step 
size dt = 0.05. Test runs with smaller dt were made to 
verify the results (see below). Total integration times 
were typically between 2 x 10 6 and 2 x 10 s units. 

Although there exist more sophisticated symplectic in- 
tegrators for Nose-Hoover baths jl5|], and although the 
standard form given in Eq.(|^) is not symplectic, the 
straightforward leap-frog was also used for them. Apart 
from simplicity and robustness, our main reason was that 
LTE should be strongly violated at the boundaries any- 
how, due to the large temperature gradients Q at the 
boundaries (see also Fig.l). If LTE would hold every- 
where, and k would depend only on T and on the sys- 
tem size, the constancy of the heat flow throughout the 
sample would imply a nearly linear temperature profile 



(since k depends smoothly on T). This is obviously not 
the case. A more direct check of LTE is obtained by esti- 
mating the normalized kurtosis (fourth cumulant) of the 
momentum distribution, 



Q, 



(<E*Lv 



(4) 



In Fig. 2 we show these cumulants for length N x = 128 
and different widths. All three curves show that indeed 
\Qi\ <C 1, verifying approximate LTE (for larger N. x we 
find even smaller Qi, as we should expect from the fact 
that temperature gradients decrease with N x ). The de- 
tailed curves depend strongly on N y , but in all cases the 
kurtosis is largest near the boundaries. Thus data ob- 
tained from the boundary regions should not be used for 
the analysis anyhow, and integration errors in these re- 
gions should not matter too much. 

Like the kurtosis, local temperatures and heat fluxes 
were indeed not calculated as functions of the spatial po- 
sitions (x, y) but for fixed particles: The temperature of 
the i-th layer is simply T, = N^ 1 (2j=i wnue the 

average heat flux through any horizontal bond connect- 
ing the layers i and i + 1 is 

Ny 

Ji = (2iV 3/ )" 1 (E(9W) + 9(i+i,j))F(ij)), (5) 
j=i 

where Fuj-\ is the force acting between particles 
and (i + l,j). We checked that indeed Ji is independent 
of i, within the statistical errors. 

Conductivities obtained by dividing the flux by the 
globally averaged temperature gradient, 



^global — 



(6) 



are shown in Fig. 3, while conductivities obtained by us- 
ing the gradient averaged only over the central half of the 
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FIG. 3: Conductivities defined by Eq.(^) for systems with 
fixed lengths N x = 32, 64, 128, . . . 8192 against N y (for N x = 
8192 there is one single point at N y = 1). Statistical and 
integration errors are less than 2 percent. 
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FIG. 4: Conductivities defined by Eq.([7|) from the same runs 
as the data in Fig. 3. 
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FIG. 5: Same data as in Fig. 4, but plotted against N x and 
with data for the same N y connected by lines. 
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FIG. 6: Conductivities defined by Eq.(]7|) for fixed aspect ra- 
tios. 



lattice, 
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N x (dT/dx) 



(7) 



with (dT/dx) centcl - = 2(T 3N:c / 4 - T Nx / 4 )/N x , are shown 
in Figs. 4 and 5. There are small but significant differ- 
ences between the two definitions of n. In general, the 
plots for K C enter show slightly more structures. We argue 
that they are more relevant, since they are less affected 
by boundary effects. We show the data for /c g i bai also, 
because they are much less noisy. 

From Figs. 3 and 4 we see that K saturates for large 
Ny. The saturation values increase roughly power- like 
with N x . For intermediate N y there are very shallow 
minima, at (N y ) nlin ~ 7V°- 35 . A rough fit to the data for 
N y > 2 is obtained with a scaling ansatz 



k = NZ'd>(N y /NP) 



a 1 = 0.26, /3 = 0.35 



(8) 



with 4>{x) « const for x 1 and <fi(x) ~ x^ a for 
x — > 0. Thus the cross-over happens at r* ~ N*~P = 
N®- 65 . As seen from Fig. 5, both the 1-d behaviour 
(r = oo) and the behaviour at r <C r* are rougly power 
like, but with different powers. In between, k increases 
even faster with N x (for fixed N y ). 

For fixed aspect ratios this leads to a decrease of k with 
N x for small N x (r > r* (N x )), and to a power increase for 
large N x (see Fig. 6). Indeed one sees definite deviations 
from a pure power law, even for large N x . This means 
that k increases asymptotically less fast than suggested 
by a simple least square fit to the present data (which 
would give the exponent a' = 0.26 used in Eq.(||)). Thus 
the true asymptotic exponent a(2d) for 2-d lattices is less 
than a'. Any detailed extrapolation is of course uncer- 
tain, but our best estimate is a(2d) = 0.22 ± 0.03. 

In order to see more precisely the behaviour of 1-d 
chains (where we have much higher statistics than for 
N y > 1), and to see more clearly the errors made by the 
finite step size dt, we show k/N x - 37 in Fig. 7. Indeed the 
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FIG. 7: The data for N y = 1, divided by iV° 37 . The upper 
two curves are for ft cen ter, the lower for K g i b a i- 

uncertainty in the definition of k is much more important 
than the integration error. Although the latter is much 
larger than the statistical error (since we obtained very 
high statistics for N y = 1), it has hardly any effect on 
the scaling of k. While the curves for Kgiobai seem to 
flatten for large N x , those for K CC ntcr seem to turn upward 
again for very large N x , after having bent downward for 
intermediate N x . Essentially the same behaviour was 
seen for the 1-d gas with alternating masses in || . Again 
any extrapolation is rather uncertain, our best estimate 
is a(ld) = 0.37 ±0.01. 

In summary, we have studied the heat conduction on 
rectangular 2-d FPU lattices with periodic lateral bound- 
ary conditions, including lattices with zero width (i.e. lin- 
ear chains). We argued that the divergence of the con- 
ductivity with system size should be universal, so that 
our results should apply also, among others, to carbon 



nanotubes. One striking result is that we do not see the 
logarithmic increase of K found in || . At least for system 
sizes up to five times the size of those studied in (§] , the 
increase is power-like, with an exponent larger than 0.2. 
This increase flattens somewhat for even larger lattices, 
but our best asymptotic estimate is still around 0.2. For 
d = 1 our estimate of the exponent a is slightly smaller 
than previous ones, but compatible with them. 

Our other striking result is that the cross-over from 
1-d to 2-d behaviour happens at an aspect ratio which 
diverges with system size. Thus it might be difficult to 
see the true 1-d behaviour in carbon nanotubes. But 
this does not mean that k should not increase with their 
length, quite to the contrary: The increase at presently 
achievable tube lengths might even be faster than both 
the asymptotic 1-d and 2-d power laws, in view of the 
special shape of the cross-over. Since k is decreasing 
with the width L± for large aspect ratios, two tubes with 
small radius conduct better than one tube with twice the 
radius. Thus, if one wants to maximize the heat flux for 
a given amount of carbon and a given length of the tubes, 
one should make the tubes as narrow as possible. 

Finally we would like to point out that the effect we 
discussed in this paper is closely related to the Casimir 
force. The latter deals with the change in energy of a 
cavity due to elimination of vacuum modes by the cavity 
walls. The enhanced conductivity of narrow strips com- 
pared to the bulk is due to the elimination of transverse 
phonons (i.e. phonons with p y ^ 0), the interactions 
with which otherwise would limit the free path length of 
soft longitudinal phonons which contribute to the energy 
transport. 

We are indebted to Roberto Livi and Antonio Politi 
for very helpful correspondence, and to Walter Nadler 
for numerous discussions. 
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